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ABSTRACT 



Context. Measuring the mass distribution of infrared dark clouds (IRDCs) over the wide dynamic range of their column densities is a 
fundamental obstacle in determining the initial conditions of high-mass star formation and star cluster formation. 
Aims. We present a new technique to derive high-dynamic-range, arcsecond-scale resolution column density data for IRDCs and 
demonstrate the potential of such data in measuring the density variance - sonic Mach number relation in molecular clouds. 
Methods. We combine near-infrared data from the UKIDSS/Galactic Plane Survey with mid-infrared data from the Spitzer/GLIMPSE 
survey to derive dust extinction maps for a sample of ten IRDCs. We then examine the linewidths of the IRDCs using 13 CO line 
emission data from the FCRAO/Galactic Ring Surrey and derive a column density - sonic Mach number relation for them. For 
comparison, we also examine the relation in a sample of nearby molecular clouds. 

Results. The presented column density mapping technique provides a very capable, temperature independent tool for mapping IRDCs 
over the column density range equivalent to Ay - 1 - 100 mag at a resolution of 2". Using the data provided by the technique, we 
present the first direct measurement of the relationship between the column density dispersion, (Tn/(n), and sonic Mach number, 
M s , in molecular clouds. We detect correlation between the variables with about 3-cr confidence. We derive the relation cr N ^ N ^ « 
(0.047 ± 0.016)A4 S , which is suggestive of the correlation coefficient between the volume density and sonic Mach number, <x p /(p) ~ 
(0.20^ 22)A4 S , in which the quoted uncertainties indicate the 3-cr range. When coupled with the results of recent numerical works, 
the existence of the correlation supports the picture of weak correlation between the magnetic field strength and density in molecular 
clouds (i.e., B oc p 5 ). While our results remain suggestive because of the small number of clouds in our demonstration sample, the 
analysis can be improved by extending the study to a larger number of clouds. 

Key words. ISM: clouds - dust, extinction - ISM: structure - Stars: formation - Infrared: ISM 
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1. Introduction 

Measuring the structure and dynamics of the dense, molecu- 
lar interstellar medium (ISM) is of crucial importance for our 
understanding of star formation from low- to high-mass stars 
and star clusters. There are a number of outstanding open ques- 
tions in star formation that come down to the coupling between 
star formation and the fundamental physical processes acting in 
molecular clouds. What determines the star-formation rates and 
efficiencies in molecular clouds and why are those parameters 
generally low? What is the relative importance of gravity, tur- 
bulence, and magnetic fields during molecular cloud evolution 
and does this vary for low- and high-mass star formation? These 
questions can be approached by studying the physical character- 
istics of molecular clouds at different evolutionary stages. Of es- 
pecial importance are studies of clouds that have not yet started 
to form stars or are in a very early stage of doing so. This is be- 
cause such clouds still bear the imprints of the processes that ini- 
tially formed them and are responsible for driving them toward 
star formation. For reviews of these topics, we refer to McKee & 
|Ostriker| ( |2007l l. 

Our knowledge on the maternal structure of clouds poten- 
tially forming high-mass stars (M > 8 M Q ) and star-clusters is 



relatively scarce. This is, in part, because of the rapid destruc- 
tion of the maternal cloud structure by stellar winds, outflows, 
ionization, and radiation pressure from the first high-mass stars 
to form in the cloud. Once the presence of high-mass stars in 
a cloud can be observationally confirmed, the maternal cloud 
structure has already been heavily affected by them. In part, the 
difficulties also originate from trivial observational restrictions: 
most sites that are thought to be favorable for future high-mass 
star and star cluster formation are located at distances of D > 1 .5 
kpc, and they are concentrated in the Galactic plane. This makes 
attaining sensitive, high-spatial resolution data on their structure 
a very challenging task. 



Only relatively recently, a class of dense molecular cloud 
clumps now commonly referred to as infrared dark clouds 
(IRDCs, hereafter) was discovered ( Perault et al.|1996{ Egan et 



|al.|1998[). These objects are currently regardec 



dates of future high-mass star formation (e.g., 



Garay et al. 2004 Beufher et al. 



2005 , Rathborne et al 



as the best candi- 



Carey et al. 



2000 



2005 
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|2006 1, although how commonly they might form massive stars 
remains still unclear (e.g., Kauffm ann & Pillai|20 1 0) > . In the most 
recent years, the study of IRDCs has seen the coming of large- 
area Galactic plane surveys that examine the Galactic-scale dis- 
tribution and properties of IRDCs through comprehensive, sta- 
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tistical samples (e.g., |Simon et al.||2006| |Jackson et al.| [2008 ; 



Peretto & Fuller|2010|]Aguirre et al.|20TT] >. 

One persistent obstacle in defining the physical conditions 
in the maternal environment of IRDCs is our inability to mea- 
sure, in a uniformly defined manner, the mass distribution over 
the wide dynamic range of (column-/mass surface-) densities 
present in IRDCs. Various observational techniques are com- 
monly used in probing the mass distributions: CO emission line 
observations (e.g., Roman-Duval et al. 2010|l, thermal d ust emis- 



sion at (sub)mm wavelengths (e.g., Schuller et al. 2009 Molinari 
|et al.|2010|, dust extinction mapping using background stars in 



near-infrared (NIR, e.g., Kainulai nen et al.|201 la| l, and dust ex 
tinction mapping in mid-infrared (MIR) using either point source 
data (e.g., Rygl et al. 2010 > or surface brightness data (e.g ., 
Butler & Tan||2009| |Ragan et al.|[2009l |Vasyunina et~aT|2 009>. 
All these mass-tracing techniques have their specific properties 
that restrict the parameter space (usually, column density range) 
within which they can provide accurate mass measurements. For 
example, CO line emission is known to be linearly correlated 
with column density only in a very narrow column density^] 
range (Ay « 3 — 6 mag) (e.g., Goodman et al.|2009"j ). Dust emis- 
sion measurements depend on the temperature of the dust grains 
that can vary by a factor of 



2 in the clouds (e.g., Bontemps 
[eTaTl[20T0l |Stutz et al.||2"0T0) |Arzoumanian et al.|[20TT] rand 
of their emissivity index. Dust extinction measurements in the 
near-infrared are limited to relatively low column densities, typ- 
ically to Ay * 1 - 25 mag in (e.g., Lom bardi & Alves]|2001 



Kainulainen et al. 2011a]). For MIR surface brightness extinction 
mapping, higher column densities can be probed up to Ay * 100 
mag, but, as discussed below, systematic uncertainties limit the 
accuracy below Ay » 10 mag. As a result, it is not straightfor- 
ward to relate the mass distribution at small, dense, star-forming 
scales to the surrounding larger-scale, more diffuse environment. 

In our previous works, we have developed techniques to 
probe the mass distributions of IRDCs by measuring the NIR 
( [Kainulainen et al.|201 la) and MIR ( |Butler & Tan|2009[ [20121 
BT09 and BT12, hereafter) dust extinction through them. While 
the NIR technique has proven to be relatively sensitive, reach- 
ing column densities of Ay m 25 mag at a resolution of ~ 30", 
the MIR data probe reliably column densities between Ay =s 
10 - 100 mag at a resolution of ~ 2". Both of the techniques 
are essentially based on a model for dust optical properties in 
infrared which are relatively well-known, and the dynamical 
ranges probed by them clearly overlap. This provides a good 
basis for a conception of combining the two techniques with 
aim of achieving both extended dynamical range and high spa- 
tial resolution. This goal is further motivated by the fact that 
the MIR technique is known to have caveats affecting especially 
low-column densities. The technique effectively filters out the 
large-scale, diffuse [Ay < 5 - 10 mag) extinction component, 
because of the estimation of the background intensity (for de- 
tails, see |Butler &"Taiil[2009l |20T2] BT09 and BT12 hereafter, 
respectively). The details of the adopted filtering scheme can 
significantly affect the performance of the technique (e.g., BT09 
Wilcock et al.|2012) . However, the diffuse dust component that 
is filtered away is well recovered by the NIR data, and it would 
seem feasible to use the NIR data to "re-calibrate" the MIR data, 



1 We note that column density (and mass surface density) are gener- 
ally related to dust extinction, Ay, through a linear relationship. To make 
comparisons with relevant earlier works easier, we use in this paper 
equivalent Ay as the main unit for column density. We give the detailed 
conversions to column density and mass surface density in Sections [2.2f 
23 



thus alleviating the uncertainty of the MIR technique arising 
from the large-scale filtering needed to estimate the MIR back- 
ground intensity. 

We develop in this paper a new technique for probing the 
mass distributions of IRDCs that takes the advantage of both 
NIR point source data and MIR surface brightness data. The 
technique is based on data that are readily available for a 
large portion of the Galactic plane, i.e., NIR data from the 
UKIDSS/Galactic Plane Survey ( [Lucas et al.|2008) and MIR data 
from the Spitzer/GLIMPSE survey ( [Benjamin et al.pOOl) . We 
will demonstrate the use of the technique by deriving column 
density data for ten IRDCs and their surroundings. 

We also demonstrate in this paper the scientific potential 
of high-dynamic -range column density data in examining the 
probability density functions (PDFs) of column densities in 
the clouds. The column density PDFs can be used to study 
the roles of gravity and turbulence in molecular clouds (e.g., 
Froebrich et al.||2007[ |G oodman ^t al.||2009] |Kainulain en et| 
|al. |2009 [ ]Froebrich & Rowle s*p010| l, generally building on the 
fact that simulations of supersonic turbulent media predict a 



log-normal shape for the density PDF (e.g., Vazquez-Semadeni 
1994; Pa doan et aL||1997a[ |Ballesteros-Paredes et al.||2011[ 
|Kritsuketal.|2011| i 



p(x) = 



0"ln 



exp 



I- 



(Inx-p) 2 

2cr2 
Inx 



(i) 



where x = p/po is the mean-normalized volume density and 
p and <xi nx are the mean and standard deviation in logarithmic 
units. Deviations from the log-normal shape are predicted for 
strongly gravitating systems, inducing higher relative fraction of 
high-column density material (e.g., Klessen 2000} [Federrath et 
|aL"1|200 8a). In K09 we presented an analysis of the PDFs of a 
sample of nearby molecular clouds and concluded that the PDFs 
of star-forming clouds are consistent with a log-normal (or nor- 
mal) shaped function at lower column densities (Ay < 2-5 mag) 
and a power-law like extension at higher columns (Ay > 2 - 5 
mag). In contrast, the PDFs of non-star-forming clouds were 
usually very well described with log-normals. This led us to sug- 
gest that quiescent clouds are primarily shaped by turbulence, 
while the effect of self-gravity is evident in the PDFs of star- 
forming clouds. However, in Kainulai nen et al.| ( |2011b| l we an- 
alyzed the data further and suggested that the PDF tail could 
be related to clumps that are not self-gravitating but rather sup- 
ported significantly by external pressure from the surrounding, 
gravitating cloud material. 

In this paper, we also use the derived column density data 
together with CO line emission data to examine the relation- 
ship between density fluctuations, measured by the column den- 
sity dispersion cr N i( N ), and turbulence energy, measured by the 
sonic Mach number Ms, in the clouds. These variables are ex- 
pected to be correlated in turbulent media and a correlation co- 
efficient between them is assumed by numerous analytic models 
of sta r formation (e.g., [Padoan et al.|1997a[ [P adoan & Nordlund 
|2lX)2l|Krumhorz & McKee|2005[|Hennebelle & Chabrier|2008t ' 
|Elmegreen|2011[ ). We measure the relationship between M s and 
°~n/(n) m the IRDCs and also in a larger sample that is sup- 
plemented by including a set of nearby molecular clouds we 
have previously derived column density data for ( Kainu lainen | 
|et al.|[2009} hereafter K09). Our results demonstrate that high- 
dynamic-range data can, and are in fact required to, probe the 
relation. 



2 



J. Kainulainen et al.: High-dynamic-range column density data for infrared dark clouds 



2. Data 

2.1. The IRDC sample 

We use as a sample the ten IRDCs selected by BT09. The co- 
ordinates of the clouds are listed in Table Q] The clouds were 
originally chosen from the sample of Rathbome et al. (2006) 
as massive, relatively nearby clouds with high contrast against 
the Galactic background, being thereby well-suited for the 8 
pm extinction mapping. The clouds are distributed between the 
Galactic longitudes of / = [18°, 53°] and kinematic distances of 
D w 1 .8 - 5.7 kpc (derived by Simon et al. 2006 ). The most mas- 



sive cloud (D) has a MIR-derived mass of M M ir = 5.34 x 10 4 
M Q . The masses are listed in Table [Tj The masses were derived 
by BT12 and they comprise the material inside the ellipsoids de- 
fined by Simo n et"aL"1 ( |2006] l, which is likely to be ~ 10% of the 
total mass of the clouds estimated over larger scales using NIR 
extinction mapping (Kainulainen et al. 2011a). Thus, the most 
massive clouds in the sample are somewhat more massive than 
the nearby Orion A molecular cloud (M * 10 5 M Q , e.g., K09). 
All clouds contain mm-cores at 11" resolution ( jRathborne et al.| 
|2006| ), most of which also appear as single cores at the 2" reso- 
lution of the MIR dust extinction data (BT12). The masses of the 
cores are typically some tens of solar masses, but depending on 
how much of the surrounding envelope is included, the masses 
can also be higher by a factor of a few (BT12). While it is cur- 
rently not well-constrained how massive a core is needed for 
high-mass star formation, the cores in principle contain enough 
material to form such stars. 



2.2. NIR dust extinction data 

We use in this work NIR dust extinction data as a tracer of low- 
extinction material (A v < 10 mag) in the IRDCs. The data were 
derived using the NIR color-excess mapping scheme presented 
by [Kainulaine n et al.|(|2011a|), which is bas ed on the so-called 
NICER technique ( |Lombardi & Alves|2001| ), but refined to deal 
with the diffuse extinction component in the Galactic plane and 
with the large number of foreground sources towards IRDCs. For 
the detailed description of the technique, we refer to Kainulainen 
et al.| (|2011a| ) and |Lombardi & Alves| ( |2001| ). |Kainulainen~et 



al. ( 201 la| see their Appendix A) presented the NIR extinction 
maps for the clouds A, B, G, H, and I. The maps for the remain- 
ing clouds were computed for this work. 

The NIR extinction mapping technique was used in conjunc- 
tion with the JHKs band data from the UKIDSS/Galactic Plane 
Survey/Data Release 6 Plus. The UKIDSS project is defined in 
|Lawrence et~aT1 ( |2007]). UKIDSS use s the UKIRT Wide Field 
Camera (WFCAM; Casali et al. 2007 ) The photometric system 
is describ ed in |Hewett et al.| (|2006| ), and the calibration is de- 
scribed in Hodgkin et al. ( 2009 ). The science archive is described 
in Hambly et al.| ( |2008[ ). The derived NIR extinction maps have 
the sensitivity (statistical uncertainty) of cr{Ay) » 0.4 - 0.7 mag, 
depending somewhat on the Galactic coordinates and the dis- 
tance of the clouds. The zero-point of the data is rather uncertain 
because of the uncertainties in the foreground source removal 
and in estimating the extinction because of the diffuse dust in 
the Galactic plane; Kainulainen et al. ( 201 la| ) give the value of 
o- (Av) ~ 1 mag for the zero-point uncertainty. The spatial res- 
olution of the data is FWHM = 30" and at that resolution they 
typically span extinction values between A v ~ 1 - 30 mag. We 
note that in derivation of the NIR extinction maps a wavelength 
dependency of optical depth, i.e., the reddening law, is assumed. 



Kainulainen et al. (|2011a| ) used the reddening law as derived by 
Cardelli et al.| ( |1989| > 



t k = 0.600t h = 0.404tj 



(2) 



It should be emphasized that the NIR technique determines 
the extinction in NIR, e.g., Ak- In order to estimate column den- 
sities, we convert the NIR extinction to visual extinction using 
the coefficient from |Cardelli et al.| ( |T989l > 



T K = 0.114TV 



(3) 



We note that we combine our results later in this paper with K09 
who used the same relative extinction law. Then, the calculated 
visual extinction can be converted to gas column density ( |Savage 
et al.|1977[|Bohlin et al.|1978||Vuong et al.|2003| ) 



Ni 



NIR 



= 2N(H 2 ) + N(H) = 1.9 x 10 21 cnT 2 (— ) 



mag' 



(4) 



Column density can further be expressed in terms of the total 
mass surface density by adding 40% (implying the H:He number 
ratio of 10:1) of helium 



2 NIR = lAmnN£ 



NIR 



(5) 



We can also express the empirical conversion from extinc- 
tion to mass surface density (i.e., Eqs. [3][5j> with the help of an 
"effective" opacity, through the definition 

T K = Kk^NIR- (6) 

Inputing Eqs.[3]|5]to this results in the effective opacity 

4 = 25.9 cm 2 g" 1 . (7) 

Corresponding effective opacities for H and J bands result di- 
rectly from Eq. [2] 



4 =43.2 cm 2 g-\ 
k) =64.1 cm 2 g _1 . 

2.3. MIR dust extinction data 



(8) 
(9) 



We use the MIR dust extinction data derived using the 
Spitzer/GLIMPSE survey 8 pm images by BT12 (see also BT09) 
as another column-density tracer for the IRDC complexes. The 
MIR mapping technique is based on measuring the decrease of 
surface brightness towards the IRDCs, which can be interpreted 
in terms of 8 pm optical depth. For the details of the mapping 
technique we refer to BT12. The derived optical depths are used 
to estimate mass surface densities, Sg, with an equation analo- 
gous to Eq.|6]and adopting a dust opacity at 8 fim 



(10) 



BT09 and BT12 chose the value k% = 7.5 cm 2 g -1 , which was 
taken as an approximate mean of five prospective dust models 
from [ Ossenkop f & Henning|(|1994|) and |Weingartner & Draine| 



(2001 1, but closest to the Ossenkopf & Henning ( 1994[ ) moder- 
ately coagulated thin ice mantle dust model (see BT09 for de- 
tails). BT09 also used weighting by the MIR background spec- 
trum and the filter response function in calculating the opac- 
ity values. Figure [T] illustrates the opacities k& and ft| tog ether 
with the opacity laws fr om |Weing artner & Draine] ( |2001| ) and 
Ossenkopf & Henning ( 1994J). The figure shows that both in 
NIR and MIR the adopted opacities represent an average of the 
shown opacity laws. It is not possible to estimate the accuracy 
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Table 1. Properties of the IRDCs. 
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Nami ng according to BT 09 



From 



Simon et al. 



(2006). 



c Saturation limit, i.e., the maximum observable mass surface density, of the mid-infrared data (cf., BT12). 

d The letters G(aussian) and F(ixed range) refer to the two approaches used in estimating the velocity dispersions in the clouds (see Section 
2.4 1. The abbreviations box and ell refer to the two approaches used in defining the cloud area (see Secti on | 3. 2b. 
e M com bined can be somewhat underestimated because of incomplete coverage of the NIR data, see Fig.|B.4 

f The distance for the cloud F has been recently derived from parallax measurements ( [Kurayama et al.| 201 1 1 and we adopt this distance 
(however, see also - 



Foster et al. 



2012} . For comparison, the GRS distance for the cloud is D = 3.7 kpc). 



The mass given by BT12 was scaled to correspond the adopted distance. 



0H94, n - 10° c 

0H94, n = initio 

WD01, R v = 3.1 

WD01, R = 5.5 

WD01, R = 5.5 Cost 



E 

" 10 




Fig. 1. Opacity laws from Ossenkopf & Henning ( 1994 1 (black 
dashed line — uncoagulated thin ice mantle dust model; black 
solid line — thin ice mantle model af ter coagulation for 10 5 
y at density of «(H) = 10 6 cm 4 ) and |Weingartner & Draine 
(f200TJ (red lines, models for Ry = 3.1, Ry = 5.5 and Ry = 5.5 
"case B" shown with solid, dashed, and dotted lines). The figure 
also shows the adopted 8 //m opacity (diamond) and the effective 
opacities corresponding to the empirical conversion of NIR ex- 
tinction to mass surface density (symbols detailed in the panel), 
see Section |Z2 



of the adopted 8 //m opacity, but the standard deviation of the 
opacities of the five models is 22%. Also, it is not straightfor- 
ward to estimate the uncertainty due to variations in the MIR 
background spectrum, but if we re-calculate the mean opacity 
using the case of flat background spectrum, a value Kg = 9.0 cm 2 
g _1 follows, indicating 20% larger value than the mean with the 
chosen background spectrum. While the flat background spec- 
trum is a rather extreme case, it seems reasonable to use the total 
uncertainty of ~30% for the absolute value of the MIR opacity. 
However, one should note that unless average dust grain prop- 
erties significantly alter from cloud-to-cloud, the relative opac- 
ity (or extinction/mass surface density) measurements between 
clouds is significantly more robust than this value. 



The Spitzer 8 /mi images have the resolution of FWHM - 
2" (1'.'2 pixel size) and that is also the resolution of the MIR ex- 
tinction data. At this resolution, the technique probes the mass 
surface density values up to £ < 0.4-0.5 g cirT 2 , or equivalently, 
Ay < 92 - 115 mag depending on the cloud. These upper lim- 
its, or "saturation" values, are estimated empirically as described 
in BT12 and listed in Table [T] The sensitivity (the absolute un- 
certainty) of the technique is not straightforward to determine, 
because the main sources of uncertainty are the variations of the 
background and foreground surface brightness. BT09 discusses 
that the technique becomes unreliable at £ < 0.013 g cirT 2 , or 
Ay < 3 mag. However, at column densities lower than 2 < 0.04 
g cirT 2 , or Ay < 9 mag, the MIR extinction maps show sys- 
tematic variations that are not present in the NIR data, likely re- 
sulting from variations in the MIR background, and hence, from 
the zero-point uncertainty of the MIR data. The MIR mapping 
technique is also affected by bright point sources and extended 
emission (nebulae) in the 8 /mi data. These were masked out 
during the mapping procedure as detailed in BT09, i.e., we do 
not estimate extinction in the regions that suffer from significant 
nebulosity. We also note that the GLIMPSE data contains some 
stripe-like artifacts resulting from non-uniform calibration. 

The properties above reflect the complimentary of the NIR 
and MIR dust extinction techniques. While the NIR data are 
more sensitive and less affected by systematic uncertainties at 
Ay < 10 mag, the MIR data can probe significantly higher col- 
umn densities. There is also an overlap range between the dy- 
namic ranges of the techniques at about Ay « 5 - 15 mag. 



2.4. CO line emission data 

We use the publicly available 13 CO (7=1-0) line emission data 
from the FCRAO/Galactic Ring Survey (Jack son et al.| |2006 
hereafter GRS) to examine the velocity structure of the IRDC 
complexes. The GRS data have a spatial resolution of FWHM - 
46", spectral resolution of 0.2 km s _1 , and typical sensitivity of 
cr(r A , ms ) = 0.4 K. 



4 



J. Kainulainen et al.: High-dynamic-range column density data for infrared dark clouds 



Since the IRDCs are located in the Galactic plane, the lines 
of sight toward them can show several emission components. We 
first determined which component corresponds to the IRDC by 
comparing their morphologies to the extinction data. This was 
done subjectively by eye, but we note that in almost all cases 
the choice was trivial; The IRDCs were originally chosen to be 
in somewhat isolated positions. In nine out of ten IRDCs the 
strongest 13 CO emission peak matched well the morphology of 
the IRDC as seen by dust extinction. Toward the Cloud E th e 



A.l 



velocity structure of 13 CO was more complicated (see Fig 
Because the dust extinction techniques can probe the total co 
umn density up to D < 8 kpc, it appears likely that the extinc- 
tion data for cloud E is affected by several, physically unrelated 
clouds. Therefore, we excluded the cloud E from the analyses 
performed in Section[4] 

We used the 13 CO line emission to estimate the velocity dis- 
persions in the clouds. The estimation was performed inside 
two regions for each cloud: an ellipsoid as defined by |Simon| 



et al. 



( 2006| l, and a box-like region that (often, but not always) 
contains a closed Ay = 7 mag contour. We describe these re- 
gions more closely in Section |3.2| and show the average spec- 
tra in Appendix [A] For both regions, we estimated the one- 
dimensional velocity dispersion, <x v , using the following two 
approaches. The first approach was based on fitting a single 
Gaussian to the mean spectrum over the velocity interval of 
the component. The velocity interval was defined using an av- 
erage spectrum over the region. We calculated the first deriva- 
tive of the average spectra, and determined the first zero-points 
of the derivative from both left and right side of the velocity of 
the maximum intensity. These zero points defined the interval 
limits. We also required the interval limit to be at most 15 km 
s _1 offset from velocity of the peak. The intervals are indicated 



A.l and A. 2 with vertical dotted lines. Then we fitted 



in Figs. 

a Gaussian to the spectrum within this interval. To get a han- 
dle on how much the velocity dispersion depends on the cho- 
sen velocity interval, we also fitted a Gaussian to the interval 
v = [vpeak - l-5<x v ,v P eak + l-5o" v ], where cr v is the dispersion 
from the first Gaussian fit. These two velocity dispersions (and 
the corresponding fits) are also shown in Figs. A.l and A. 2 and 
we give in Table[T]their mean value. The mean value is also used 
in all further analyses. 

The second approach for calculating the velocity dispersions 
was calculating the standard deviation of the spectrum in the 
fixed velocity interval whose definition is described above. As 
the mean in this calculation, we used the intensity-weighted 
mean, which represents the center-of-mass in the case where in- 
tensity is linearly proportional to the column density along the 
line of sight. The linear proportionality assumes that the 13 CO 
emission is optically thin and does not suffer from depletion (see 
|Hernandez & Tan]|2011| [Hernandez e t al.||20fT for a study of 
these effects in the clouds F and H). Only channels with intensity 
higher than three times the root-mean-square noise level were in- 
cluded in the calculation. The resulting velocity dispersions are 
given in Table [T] 



In addition to the GRS data, we use in Section 



5.2.2 12 CO 



data from the COMPLETE survey ( |Ridge et al.||20Q6) > and the 



survey of Dame et al. ( 200 \\ to examine the velocity dispersions 
of nearby molecular clouds. The COMPLETE data was avail- 
able for Ophiuchus and Perseus. The Dame et al. (2001) data 
was adopted for the clouds in the sample of K09 (except for the 
Lupus clouds and LDN134 that are not covered by Dame et al.). 
The CO data for these nearby clouds do not suffer from the line- 
of-sight confusion (they are outside the Galactic plane) and the 
linewidths toward them were computed from a Gaussian fit to 
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Fig. 2. Pixel-to-pixel comparison of the NIR and MIR column 
density data for the cloud B. The dashed line shows the one-to- 
one correlation. The data show correlation below Ay lR < 10-15 
mag, above which the NIR data starts to underestimate extinc- 
tion compared to the MIR data. The offset of ~ 5 mag in the 
range A^ IR < 10 mag originates from the spatial filtering of the 
MIR technique. 



the mean sp ectrum. In addition, we adopt the mean linewidth for 
Taurus from|Goldsmifh et al.|(|2008|, i.e., cr( I3 CO) = 1.0 km s" 1 . 



3. Methods 

3. 1 . Combination of NIR and MIR data 

We combine the NIR and MIR dust extinction maps with the aim 
of achieving high-dynamic-range extinction data. We first exam- 
ined the correlation between the data sets by convolving the MIR 
data (FWHM = 2") to the significantly lower resolution of the 
NIR data (FWHM = 30") and making pixel-to-pixel compar- 
isons. As an example, Fig.[2]shows the comparison for the cloud 
B. Despite the rather high scatter, the data sets correlate at col- 
umn densities below Ay < 10-15 mag. In this range, the column 
densities measured by the MIR data are systematically about ~4 
mag lower than the column densities measured from NIR data. 
The offset is present in all clouds and varies between Ay ~ 2 - 8 
mag. The offset is most likely a result of the spatial filtering of 
the MIR mapping technique which effectively subtracts the dif- 
fuse cloud component from the MIR maps. At column densities 
higher than A v > 10-15 mag the NIR data begins to underes- 
timate the column density, which results from the fact that the 
NIR data does not efficiently sample high (Ay > 30 mag) ex- 
tinction lines-of-sight (none, or very few, background stars are 
visible through higher-extinction clumps/cores). 

Since the NIR and MIR data showed correlation in the col- 
umn density range Ay < 10-15 mag, we conclude that it is 
reasonable to assume that the difference between the MIR and 
NIR extinction data below this range is a good representation of 
the background component that the MIR technique filters out. 
With this in mind, we combined the NIR and MIR data with 
the following procedure. First, we re-binned and smoothed the 
MIR data to the resolution of NIR data (FWHM - 30", pixel 
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size 15") and subtracted it from the NIR image to generate "a 
correction image" 



A correction (3(r) = ^NIR _ A MIR( 30 »)_ 



(ID 



When built in this way, the correction image, Ay mctlon (30"), 
approximates the large-scale component filtered out by the MIR 
technique at spatial scales equal or larger than the background 
filter size used in the MIR technique (in the case of BT12, a few 
arcmin). Thus, the correction image clearly resolves the filtering 
scale of the MIR technique. Since the NIR data underestimates 
the column density above Ay > 10-15 mag, we blanked the 
pixels in the correction image in which the extinction is higher 
than Ay > 10 mag. Values for these pixels were estimated as 
a mean within an aperture of 75" (5 pixels of the NIR data). A 
Gaussian weighting function with FWHM = 30" was used to 
weight the values when calculating the mean. Then, the correc- 
tion image ^^ olrectlon (30") was re-binned to the pixel scale of the 
MIR data (l'/2) and added to the original MIR column density 
map 



\ combined _ \ NIR , \ i 

/ly — /iy + /iy 



correction 



(2"). 



(12) 



This resulted in column density maps, pf° , that take advan- 
tage of the signals both in the NIR and MIR data. The pixels that 
lack either NIR or MIR data (e.g., because of the incomplete 
coverage of the NIR data, bright MIR point sources or extended 
MIR nebulosity) were flagged to have no data in the maps. As a 
result, extinction is not determined for every pixel but only for 
regions in which we consider both the NIR and MIR techniques 
to be unaffected by any artifacts. We describe the results further 
in Section 



4.1 



In the remainder of this paper, when we refer 
to dust extinction/column density/mass surface density maps, or 
data, we refer to these combined maps. 

We finally note that the combination procedure is not sensi- 
tive to the absolute values of the adopted dust opacities (see Fig. 
[TJ, only to the ratio of dust opacities in NIR and MIR. Only the 
scaling of the observed extinctions to the units of mass surface 
density requires knowledge of the absolute dust opacity. In our 
approach the NIR-MIR dust opacity ratio equals to 



A'8 



0.29 ±0.10. 



(13) 



As shown in Fig. [T] this ratio represents an average behavior of 
the displayed models. We quote as uncertainty of the ratio the 
standard deviation of the models used in Fig. [T] which is 30% 
for the 8 pm (see Section |2~3] l band and 15% for the K band. 
Combining these values quadratically results in uncertainty of 
36%. However, we again point out that this measure of uncer- 
tainty reflects a selection of dust models. If the dust grain prop- 
erties are relatively similar from cloud-to-cloud, the ratio in con- 
siderably less uncertain. Currently, there is no detailed infor- 
mation available on this ratio in the particular environment of 
IRDCs, and hence, aligning the relative reddening law more ac- 
curately based on actual measurements is not possible. 

3.2. Definition of the cloud area 

It is not straightforward to define an object such as "a cloud" in 
the relatively crowded fields toward the Galactic plane in which 
the IRDCs reside. In Sections [5. 2.1 115. 2!2l we examine statistics 
derived from the column density PDFs of the IRDCs. Despite the 
seemingly straightforward definition of the PDF, it is not trivial 
to measure its shape from observational data in a robust way. An 



ideal target would be a totally isolated cloud in which the lowest 
column density included in the PDF forms a closed contour in 
the column density map. Such a closed contour can be regarded 
as the completeness limit of the PDF; The shape of the PDF be- 
low the level of the lowest closed contour can clearly suffer from 
incompleteness. For example, in K09 we used as a definition of 
a cloud the level Ay = 1 mag, which in most cases was below 
the mode of the PDF. In the case of the IRDCs that are located 
in the Galactic plane the confusion from physically unrelated 
clouds that are nearby in projection makes it more difficult to 
define the region which should be included in the PDF. 

With this caveat in mind, we adopted two definitions of the 
cloud area for our analysis. The first definition is to use the ellip- 
soids defined by Simon et al. (2006). We chose from the catalog 
the ellipsoids with largest area within the mapped area, which in 
practice roughly encircles the high-column density (Ay > 10-15 
mag) regions in the IRDCs. The ellipses are shown in the figures 
of Appendix [5] The second definition of the cloud area was con- 
structed with the aim in being as complete as possible down to 
relatively low column densities. In the second definition, we de- 
fined a rectangular area that surrounds the Simon et al. (2006) 
ellipsoid, and makes an effort to completely surround a closed 
Ay - 1 mag contour. However because of blending of neigh- 
boring clouds in two-dimensional data, this contour is not nec- 
essarily closed in all clouds. In cases where the blending was 
severe, we used the velocity structure as traced by the 13 CO data 
to help disentangling physically unrelated components. The final 
rectangular regions chosen are shown in the maps in Appendix 
IB] We conclude that the data within these areas should be close 
to complete above Ay > 7 mag. However, we also note that in 
some, more isolated clouds slightly off the Galactic plane, such 
as clouds C and J, the lowest closed contour can be as low as 
Ay x 2 mag. For brevity, we refer to this second cloud definition 
as a "Ay = 7 mag box" in the remainder of the paper. 



4. Results 

4. 1 . High-dynamic range column density maps of IRDCs 



We used the method described in Section 13711 to derive column 
density maps for the ten IRDC complexes of our sample. We 
present the maps in Appendix [B] (online only). The pixels in the 
maps for which there was either no MIR or NIR data available 
are blanked to black. These blanked areas can result from strong 
MIR point sources and nebulosity, or missing UKIDSS data (due 
to the incomplete coverage of the survey). 

Figure [3] shows as an example the column density map of 
the IRDC H and, for comparison, the separate NIR- and MIR- 
derived column density maps of the complex (Kainulainen et al. 
|201 la| BT12, respectively) in the same color-scale. The com- 
bined map clearly reflects the characteristics of both data sets: it 
recovers the large-scale structures in the cloud with better sen- 
sitivity than the MIR data alone, and it retains the high angu- 
lar resolution of high-column density structures of the original 
MIR-derived maps. The combined maps probe column densities 
from the sensitivity limit (cr(A v ) w 1 mag) to the saturation limit 
of the MIR data (A v , sa t - 100 mag). 

We list in Table Q] the new masses of the IRDCs calculated 
from the combined column density data. The masses from the 
combined data are on average 1 .4 times larger than the masses 
calculated from the MIR data alone. We also show in Table Q] 
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also the virial parameters of the clouds, defined as ( |Bertoldi"&| 
|McKee|1992> 



1 .0000 



5o- 2 v R 
GM ' 



(14) 



In this <t v is the velocity dispersion in km s _1 for which we use 
the dispersions derived using Gaussian fits and ellipsoidal cloud 
definition (i.e., 0"o,eii m TableJTJ. R is the radius of the ellipsoid in 
pc for which we use the expression R - y/aiea/n:, and M is the 
mass in solar masses. G = 1/232 is the gravitational constant 
in corresponding units. Five out of ten clouds have their virial 
parameters within the factor of 2 of unity. The cloud G has the 
largest virial parameter, a wu — 4.7. None of the clouds have their 
virial parameters more than factor of 2 below unity. Modulo the 
effects of surface pressure terms, internal density structure, and 
cloud elongation, these results suggest that IRDCs (at least those 
typical of this sample) tend to be self-gravitating, gravitationally 
bound objects (see also |Hernandez et al.|2012| l. 

4.2. Probability distributions of column densities 

We used the derived maps to examine the column density PDFs 
in the IRDCs. The PDFs describe the probability of observing a 
column density value Ay in the clouds. Because of the complete- 
ness issues in defining the clouds with the Sim on et al.| (2006) 
ellipsoids, we only examine the PDFs within the A v = 7 mag 



boxes (see Section 3.2 1, and analyze them only above Ay > 1 



mag. We also note that we do not use all the dynamic range 
present in the column density maps which extend typically to 
Ay - 140 - 150 mag. As discussed by BT12 the saturation of 
the mid-infrared technique starts to affect the column density 
values clearly earlier. The saturation limits, L sat , below which 
BT12 conclude the maps to be relatively unaffected by satura- 
tion, are given in Table [T] The error in column density due to 
saturation is for most clouds < 10% when S < 0.82 sat . We only 
include the column densities up to the saturation limits S sat to 
the PDFs. 

Figure [4] shows the PDFs of eight IRDCs of our sample. The 
PDFs were normalized for the plot by first fitting a line to the 
range In Ay = 2-3 (Ay = 7.4 - 20 mag) and then setting a 
boundary condition p(\nAy = 2.0) = 1. In general, the PDFs of 
all IRDCs are very similar to each other below Ay < 30-40 mag. 
Above that, they differ from each other, with the PDFs of some 
clouds decreasing more rapidly than others (i.e., they contain a 
lower fraction of high column densities). For reference, we show 
in Fig.|4]a log-normal function with fi = 1.2, <xi nAr = 0.9 (see 
Eq. []}. The PDFs of all but two IRDCs (C and D) are within a 
factor of three of this curve throughout their dynamic range. The 
individual clouds show varying behaviors, e.g., the cloud C has 
almost a power-law like PDF until Ay « 70 mag, after which 
it drops abruptly. It is not easy to assess whether the PDFs are 
better fit with a power-law between In Ay m 2-3.5 and a drop 
afterwards, or by a log-normal. This is because of the similarity 
of the power-law and a wide log-normal function whose peak 
(mode) is not covered. We analyze the properties of the PDFs 
further in Sections [5.2.1115.2.21 

We also note as a caveat that the pixels affected by bright 
MIR sources (YSOs) and nebulae are blanked in the extinction 
data and can affect the PDFs (see Section 3.1 1. In principle, it is 



possible that especially the high-extinction portions of the PDFs 
are affected by the missing data, because high-column density 
regions could presumably harbor more active star formation than 
low-column density regions. There are bright MIR point sources, 




Fig. 4. Column density PDFs of the IRDCs in the sample. The 
PDFs are normalized so that a linear fit to the range lnA v = 
2-3 fulfills the boundary condition p(lnAy = 2.0) = 1. The 
dashed histogram shows the PDF of the Ophiuchus cloud from 
K09. The black dotted vertical line marks the upper limit of the 
K09 data (Ay = 25 mag). The black dotted curve shows, for 
reference, a log-normal function with fi = 1.2, <xi nN = 0.9 (not a 
fit). The colored vertical dotted lines show the saturation limits 
of the IRDCs. 



likely local to the cloud, present inside the Simon et al. (2006) 
ellipsoid areas of the clouds A, C, D, E, and H. However, the area 
covered by these blanked regions is only a fraction of, e.g., the 
area above Ay > 40 mag, which gives some suggestion that their 
effect should not be very large. Similarly, while almost every 
cloud has some MIR-bright nebulosity in their area, only in the 
case of Cloud I is the fraction of blanked pixels higher than a 
few percent (for the cloud I about 10% of the "7 mag box" is 
affected, see Section [3~l2"| >. Therefore, we note that some caution 
is in place when regarding the PDF of the Cloud L and we also 
exclude it from the analyses of the PDFs performed later in this 
paper. 

4.3. Fractional amount of high-column density material 

We examined the fractional amount of high-column density ma- 
terial in the IRDCs by calculating the cumulative distribution 
functions (CDFs) of their column densities. Here we define the 
CDF as the fraction of mass above the reference value 



dM'(> A v ) = 



M(> Ay) 
Mjot 



(15) 



where M(> Ay) is the mass above the threshold Ay and Mtot 
is the total mass above A v = 7 mag. The CDFs are constructed 
from the same pixels as the PDFs in Section 4.2 i.e., for the 
Ay = 7 mag box. The resulting CDFs are shown in Fig. [5] (left 
panel). The CDFs show a relatively linear shape in the log-linear 
scale, i.e., their shapes are consistent with exponential functions. 
We find that on average, the CDFs of the clouds follow the form 



dM'(> A v ) cc e 



-0.07X/W 



(16) 



This relation was obtained through a fit to the mean CDF of all 
IRDCs. However, we note that the CDF, in general, depends on 
the dynamic range and resolution covered by the observations. 
Remembering this, we compared the CDFs derived for IRDCs 
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Fig. 3. Column density maps of the H cloud. The contours in the three-panel figure are at Ay = {7,40} mag. The color scale is given 
in units of visual extinction, Ay (black numbers), and mass surface density, 2 (red numbers). Top: The map based only on NIR 
data (Kainulai nen et al.|201 la| l. Middle: The map based only on 8 um data (BT12). Bottom: The map based on combined NIR and 
MIR data. The standard deviation of column density values in a region about ~ 10' distance from the cloud where Ay < 2 mag is 
cr(Ay) = 1.4 mag. The maximum values in the map are Ay w 140 mag. The separate panel at the bottom shows a blow-up of the 
combined extinction map. The contours are at A v = {7,n x 10} mag. 
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with those of nearby molecular clouds from K09. To facilitate 
the comparison, the dynamic ranges of the IRDC data and the 
data from K09 were truncated to the range covered by both data 
sets. The lower limit for this range is set by the IRDC data, which 
is complete only above Ay > 7 mag. The upper limit is set by 
the data of K09 that probe only values up to about Ay < 30 mag. 
Therefore, we constructed the CDFs for this comparison using 
only the range Ay = 7 - 30 mag. 

Figure [5] (right panel) shows the resulting CDFs for the 
IRDCs and three nearby giant molecular clouds (Orion A, Orion 
B, and the California Cloud). These nearby clouds have similar 
masses (M * 10 5 M ), but their star-forming activities differ 
significantly. Orion A is the most active cloud with 2862 mem- 
ber sources, Orion B has 635 members, and the California Cloud 
only 279 members (adopting the member counts from |Lada et al 



2009) . Fits ofEq.[T6 
0.18 



-0.22, 



to the CDFs of these clouds result in slopes 
and -0.49 for the Orion A, Orion B, and California 
Cloud, respectively. The average CDF of the IRDCs has an ex- 
ponent of -0.12, which is even flatter than that of Orion A. The 
main source of uncertainty in the CDF of the IRDCs originates 
from the chosen absolute value of k& (see Section [3~Tj ). To first- 
order, the uncertainty in Kg reflects as a multiplier to the CDF: for 
example, underestimation of Kg by a factor of 2 results in over- 
estimation of column densities by a factor of roughly two, and 
hence, the exponent of the CDF will be under-estimated by the 
same factor. In Section l2~3l we discussed that the absolute uncer- 
tainty of k$ is ~ 30%. This translates to uncertainty of ~ 30% 
for the exponent of the CDF. Hence, the exponent -0.12 we de- 
rive from the IRDCs is roughly 2-cr higher than that of Orion A. 
We conclude that this implies that the fraction of high-column 
density material in IRDCs is at least comparable, and possibly 
higher than in nearby active star-forming clouds. This result is 
similar to what we found out in Kainulainen et al. ( 201 la| > using 
only NIR data. With the data presented in this paper we quantify 
the form of the CDF above Ay > 7 mag and establish that it on 
average has an exponential shape that continues up to Ay < 100 
mag. 



5. Discussion 

5. 1 . Column density probability distributions of IRDCs 

How do the PDFs we derive for IRDCs in this paper compare 
with earlier findings of the column density PDFs in molecular 
clouds? We show in Fig. [4] the PDFs of the IRDCs and that of 
the Ophiuchus cloud from K09, which is representative of most 
star-forming clouds. The PDF of Ophiuchus is very similar to 
those of the IRDCs, and consistent with either a power-law or 
a log-normal Ou = 1.2, cr toN = 0.9). Thus, the PDFs of IRDCs 
seem to resemble the PDFs of nearby star-forming clouds, al- 
though they are slightly shallower, i.e., they contain relatively 
more high-column density material. 

We note that in K09 we refer to the high-column density 
tails (at Ay > 2 - 5 mag) of the PDFs as "power-law" tails. 
We elucidate the nomenclature and the dynamic ranges of the 
K09 data and this work in Fig. [6] The power-law-like behavior 
of the PDF of Ophiuchus is obvious in Fig. |4] However, the PDF 
of Ophiuchus is, as well as the PDFs of the IRDCs are, consis- 
tent with the log-normal function shown in the figure. This log- 
normal function is significantly wider than those used by K09 to 
describe the PDFs of low-column density (Ay < 2 - 5 mag) ma- 
terial. Thus, the data we present in this paper suggests that also 
the high column densities in molecular clouds can be described 
by a log-normal PDF, however, one clearly wider than the log- 
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Fig. 5. Left: Cumulative distribution functions (CDFs) of the 
IRDCs, expressing the fraction of mass at column densities 
higher than the abscissa value. The functions are normalized to 
unity at Ay =7 mag. For comparison, the CDF of the Orion 
A molecular cloud from K09 is shown in the panel. The func- 
tions are normalized at their value at Ay = 7 mag. Right: The 
CDFs of the IRDCs and three nearby clouds from K09 (Orion 
A, Orion B, and the California Cloud, shown with dotted lines). 
The CDFs of the IRDCs were derived after smoothing and re- 
binning the IRDC column density data to the spatial resolution 
of the nearby clouds. Only the values in the range Ay =7-30 
mag were used in deriving the CDFs shown in the panel. 
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Fig. 6. Column density PDFs of a nearby star-forming cloud 
(Taurus, red histogram) from K09 and one IRDC of this study 
(Cloud A, blue histogram). In K09, the different parts of the 
PDFs are referred to as "log-normal" and a "tail" (or, a "power- 
law)". The ranges of these two components are indicated in the 
figure with red arrows. The red dashed line shows a log-normal 
fit to the "log-normal" part of the PDF. The blue dashed line 
that shows a log-normal fit to the IRDC data results to a clearly 
wider log-normal function than those discussed by K09. As the 
figure indicates, the log-normals we fit to IRDC data are rela- 
tively close in their shape to the "power-law tails" discussed by 
K09 in the range Ay = 7 - 25 mag. 



normal that we used in K09 to describe low column densities. 
It remains to be explored in future studies whether the column 
density PDFs of IRDCs are well-described by log-normals down 
to low column densities (A v m 1 mag), or whether they show a 
similar change in shape at around Ay = 2 - 5 mag as the nearby 
star-forming clouds do. 
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In K09 we showed that the non-star-forming clouds do 
not show similar high-column density tails as the star-forming 
clouds do (or, considering the discussion above, their PDFs 
are clearly narrower than those of the star-forming clouds). 
This could be interpreted as gravity having a significant effect 
in shaping the cloud structure as has been established in nu- 
merical studies (e.g., [Klessen 2000; Ballesteros -Paredes et al.| 
201 1) , but also possibly as confinement of non-self-gravitating 
clumps by the pressure imposed to them by their surrounding 
medium (Kainulai nen et alT||2011b| l. What is common in both 
of these frameworks is that the shape of the PDFs, and espe- 
cially the presence of a wider, log-normal (or power-law-like) 
component, is intimately linked to the onset of star formation. 
Further, Rowles & Froebrich ( 201 l| l analyzed the PDF shapes 
and embedded populations of nearby clouds and suggested that 
the PDF shape is also linked to the star-formation mode (iso- 
lated/dispersed vs. clustered) of the clouds. The PDFs of the 
IRDCs resemble more closely the PDFs of the star-forming 
clouds than non-star-forming clouds (cf., Figs. |4j [5]) in a sense 
that the fraction of high-column densities is high (regardless 
whether the PDF shape is log-normal or power-law). This alone 
suggests that substantial (low-mass) populations may be present 
in the IRDCs of our sample. 

We established in Section 14.31 that the fractional amount of 
high-column density material in IRDCs in the range Ay =7-30 
mag at least comparable to, and possibly even higher than in 
active nearby clouds. While there is a well-known correlation 
between the gas surface density and star-forming rate (SFR; 
Sch midt||l959| |Kennicutt|[l998| ), it was recently proposed by 
Lada et al.| ( |2012| l that there exists a more fundamental scaling 
between the surface density of dense gas and SFR. In particu- 
lar, Lada et al. ( |2012[ ) suggested that the surface density of SFR 
depends linearly on the mass surface density of H2 



£SFR <* /dg£(H 2 ), 



(17) 



where /dg is the fraction of dense gas in the clouds, defined as 
M(A K > 0.8 mag) 



/dg - 



M(A K > 0.1 mag)' 



(18) 



Here, Ak is the extinction in the K band (Ak = 1 mag is equiv- 
alent to A v ~ 8.8 mag, Eq. |2j. The fraction /dg was measured 
and analyzed in a number of nearby clouds by Lada et al. (2010 1, 
who derived the mean value (/dg) = 0.1 ± 0.06, with maxi- 
mum fractions of / D g ~ 0.2 in the Orion A and Corona Australis 
clouds. Because the PDFs of the IRDCs we present in this paper 
are complete only above Ay > 7 mag, we cannot directly calcu- 
late the values of / D g from our data. However, we can present an 
estimate for the /dg in the following manner. We re-examined 
the CDF data presented by K09. We fitted exponential func- 
tions to the high-column density slopes of the CDFs of all star- 
forming clouds (13 clouds) and extrapolated the clouds' masses 
above Ak > 0.1 mag using the fitted functions. This test showed 
that the cloud masses were underestimated on average by a fac- 
tor of (c) = 1 .9*q! using this procedure. The scatter in the factor 
had the 3-cr interval of c = [0.7, 5.2]. 

Since the cloud masses are systematically underestimated, 
we applied the factor (c) = 1 .9 to the masses of IRDCs derived 
by extrapolating their exponential CDFs to lower column den- 
sities. This resulted to the average dense gas mass fraction of 
/dg ~ 0.2 of our IRDC sample, which is about two times larger 
than the average fraction derived for nearby clouds by Lada et al. 
(2010). Given the uncertainty of our analysis, we find that this 
implies that the /dg in IRDCs is of similar level, and possibly 



even higher, than in nearby clouds. This, in turn, suggests that 
the SFR of the IRDCs should be at least comparable to equally 
massive nearby molecular clouds on average. We recognize that 
this prediction remains suggestive, because we could not directly 
calculate the dense gas mass fraction, but had to rely on a cor- 
rection for the missing mass below Ay < 7 mag. 

5.2. Relationship between the sonic Mach number and 
column density variance in molecular clouds 

We demonstrate in this section the potential of high-dynamic- 
range column density data to examine the relationship between 
the column density dispersion and the sonic Mach number in 
molecular clouds. This relation is interesting, because it is di- 
rectly linked to the fundamental relation between density fluctu- 
ations and turbulent energy in supersonically turbulent media: 
schematically, higher turbulent energy (higher Mach number) 
is expected to induce stronger density fluctuations, and hence, 
wider density PDF. In the following, we describe the formalism 
of this connection in more detail. 

Numerical simulations predict a log-normal form for the 
PDF of volume densities in non-gravitating, isothermal media 
(e.g. Vazquez-S emadeni 1994 Padoan et al. 1997a). The dis- 
persion of the logarithmic standard deviation cri nx (see Eq. [T| is 
linked to the turbulent energy, i.e., the (non-thermal) sonic Mach 
number, M s , via the equation (e.g., Nordlund & Padoan 1999; 
IQstriker et alTl999l ) 



G^ nx =ln(l+M^ 2 ), 



(19) 



where b is a proportionality constant. This equation translates 
to the dispersion of non-logarithmic, mean-normalized volume 
densities 



r p/(p> 



bM s . 



(20) 



Thus, the proportionality constant b concretely describes how 
strong is the coupling between the turbulent energy and density 
fluctuations. Exploring the value of b in numerical simulations 
(e.g., |Padoan et al.||1997al JKritsuk et al.||2007} |Federrafh et al. 
2008b; L emaster & Stone||2008[ |Price et al.||201 l| l and analyti- 
cally (e.g., Molina et al.|2012j has received ample attention re- 
cently, owing to its intimate connection to the analytical formu- 
lations of the dense core formation rate and the possible origin 
of the IMF (e.g.,|Padoan & Nordlund|2002|[K rumholz ~& McKee 
[2^03tpennebelle & Chabrier|2008||Elmegreen|2011| >. 

While observations of column density only probe the 2D 
density statistics, it is established that the log-normal nature 
of the volume density PDF remains in 2D to 3D transforma- 
tion (e.g., I Vazquez-Semadeni & Garcla]|2001| |Federrath et al. 
|2010[ ). Taking this statement one step further, |Brunt et aL ( 2010) 
presented a formalism to directly estimate the 3D density vari- 
ance from the measurements of 2D column density field. Their 
method is based on the possibility to estimate the 3D-2D vari- 
ance ratio 



NRN) 



CT 



(21) 



pl(p) 



using solely observations of 2D density field. In the equation 
o~n/(n) is the standard deviation of mean-normalized column 
densities. We note that Brunt et al. (2010) examined the values of 
R in a set of isothermal simulations of supersonic turbulence and 
found that in all their models R was between R = [0.03, 0.15]. If 
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we now combine Eq. [2T]with Eq.[20] we have an expression for 
b which depends only on 2D statistics 



Ms M s 



(22) 



In this equation, ci\ is the ratio between the sonic Mach number 
and column density variance, i.e., the slope potentially measur- 
able (however, see the discussion below) through a linear fit to 
the [Ms,cr N / W ] data 



(23) 



It remains to be explored in future work whether the quality of 
our combined NIR+MIR data is suitable to directly estimate the 
variance ratio R, and from therein, the 3D variance in IRDCs. 
Nevertheless even in the absence of the transformation, rela- 
tionship between the 2D density variance with the sonic Mach 
number is a proxy of the relationship in 3D. Supporting this 
approach, recently Bur khart & Lazarian| ( |2012| > showed that in 
MHD simulations the Ms - <T\ BX relation is preserved in trans- 
formation from 3D to 2D variables, with the details of the corre- 
lation depending somewhat on the simulation parameters. 

One should emphasize that interpreting observational 
(M s ,o-n/(n)) data in terms of the formalism presented above 
assumes that the measurements of the column density variance 
and the mean column density are performed over such dynamic 
range that they represent the true underlying total variance 
and mean. This is clearly not the case if the column densities 
are measured only over a relatively narrow dynamic range, as 
would be the case with, e.g., 13 CO line emission observations 
(A v * 3 - 6 mag, e.g., |Goodman et al.|200"9] >. 

With this background in mind, we present in the following 
two approaches for analyzing the Ms - ctnkn) relation. First, we 
describe a calculation in which we simply correlate the column 
density dispersions measured from the Ay data above Ay > 7 
mag with M s . This experiment can serve as a point of com- 
parison for other numerical/observational works, but we cannot 
directly connect it with the formalism presented earlier in this 
section. Second, we make an effort to estimate the mean and 
variance of the observed column density fields in a dynamic- 
range-independent manner by assuming that the underlying col- 
umn density PDF is log-normal. This assumption allows us to 
connect the relation directly to Eqs. 20p3 and to estimate how 
strongly the density fluctuations couple with the turbulent en- 
ergy, i.e., estimate the coefficient b (Eq.|22|). 

5.2.1. Model-independent approach 

We first examine the relationship between the sonic Mach num- 
ber, M s , calculated from the 13 CO line emission and the standard 
deviation of the mean-normalized column density data, ct^tVa. 
defined af] 



data 

N I (N) 



1 " 

(=1 



\(N)'> 



(24) 



where n is the number of column density values N that enter the 
calculation and (N) is the mean column density in the field. We 
examined Ms and 0"^'^ only in the case of the Ay = 7 mag 



2 We use here a notation ir™™ to emphasize that the parameter is 
calculated directly from the column density data and is not necessarily 
the same as cr N /( N ) discussed in Section 5.2 



box, because of potential completeness issues in o^ff^ with el- 
lipsoidal cloud definition. Only data above A v > 7 mag was 
considered. The Mach numbers were calculated from the one- 
dimensional velocity dispersions derived in Section |2.4| (given 
in Table [T} with the expression 



Ms 



3D 



CO 

Cs 



V34°o 



(25) 



where c s = (kT ///mn) 1 ' 2 is the isothermal speed of sound and 
the velocity dispersion of CO is assumed to represent turbulent 
velocity. In the lack of measurements of the gas kinetic temper- 
atures, we chose T = 15 K, consistent with the observations 
of |Carey et aT] ( fT998] l and |Pillai et aL] p006| , and n = 2.33 
for the calculation of the speed of sound for all IRDCs. Eq. 25 
assumes that the velocity field is isotropic and the 3D velocity 
dispersion is related to the line-of-sight velocity dispersion with 

0~3D = V3c 1D . 

Figure |7j shows the relation between M s and c^j^ in the 
IRDCs. The figure shows the data separately for the two dif- 
ferent methods to calculate Ms (see Section |2.4| >. The Mach 
numbers of the clouds span the range Ms = 8-35, while 
the column density dispersions extend over a narrower range of 
cr ty(N) ~ 0-4 - 0.8. This range is in agreement with BT09 who 
concluded that there are differences in the widths of the column 
density PDFs of these clouds (BT09 used the 'ellipsoid' area in 
determining their cloud region; see Section 3.2 1. The small num- 



ber of clouds in our sample limits the ability to establish whether 
the two variables or correlated or not. Also, the uncertainties 
in both variables are likely dominated by the choices made in 
defining the cloud region instead of statistical uncertainty. The 
Pearson's correlation coefficients for both choices of velocity 
dispersions are below 0.8 (which corresponds to the p-value of 
p w 0.02 when degrees of freedom is 6). We also fitted the data 
with the linfitex routine from the mpfitfun library (Mar kwardt 
2009 ). In the fitting, we used as uncertainties the relative uncer- 
tainty of 20% for o^fjv) (originating mostly from the zero-point 
error of the data, cr(Ay) w 1 mag) and 30% for M s (originating 
mostly from the uncertainty in temperature, which is typically 
10-20 K in IRDCs). The resulting slopes (and their errors) im- 
ply that the correlation of variables is not significant. However, 
we emphasize here that, as explained in Section 5.2 this relation 



cannot be straightforwardly linked to the theory describing the 
M s - o-pi(p) relation in turbulent media. We also note that uncer- 
tainty in opacity-law somewhat increases the uncertainty of the 
derived relationship (see discussion in Section 5.2.2 1. 



5.2.2. Approach based on a log-normal column density PDF 

We can connect the observed column density PDFs and M s val- 
ues intimately into the theory of turbulent media (Section 5.2 1 by 
interpreting the PDFs shown in Fig.|4]as the high-column den- 
sity sides of log-normal-shaped PDFs and estimating the shape 
parameters of the underlying PDFs through fits of log-normal 
functions. The log-normal function is given by Eq. [T[ and we 
use in it the variable x — N, which is the (non-normalized) col- 
umn density. With this definition, ju and cr ln N are then the mean 
and standard deviation (in logarithmic units) of the log-normal 
distribution of column densities. Note that this formulation de- 
scribes the PDF of non-normalized column densities: we do not 
know a priori the underlying mean of the column density field 
which should be used for normalization. From the fit, we deter- 
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Fig. 7. Standard deviation of the mean-normalized column densities, cr^^, as a function of the sonic Mach number, M s , in eight 
IRD CS. The Cloud E is not included in the plot, because of the potential problems in the extinction data of the cloud (see Section 
2.4 1. The Cloud I is also excluded because of nebulosity in the cloud region. The left frame shows the relation when /W s is calculated 
from Gaussian fits to the mean spectra. The right frame shows the relation when yV( s is calculated from a fixed velocity range (see 
Section 2.4 1. The Ay = 7 mag box was to define the cloud region that was included in the plots (see Section 3.2 1. 



Table 2. Parameters of the log-normal fits to the PDFs 



Cloud 


(N) [mag] 


o- N [mag] 


O" N/(N) 


A 


4.1 


5.3 


1.3 


B 


6.1 


4.4 


0.7 


C 


8.3 


13.7 


1.7 


D 


6.7 


7.5 


1.1 


E 








F 


4.3 


8.9 


2.1 


G 


5.5 


6.1 


0.9 


H 


4.9 


4.2 


0.9 


I 








J 


8.8 


6.2 


0.7 



mine the parameters (p., cr] nN ), which are used to calculate the 
mean, (N), and standard deviation, <x N , in non-logarithmic units 



(N) 



cr N = [(e' r r„N _ l) e W„N] 



1/2 



(26) 



(27) 



The standard deviation of the mean-normalized, non-logarithmic 
column densities is finally the standard deviation of non- 
normalized data divided by the mean m 



cr N 



(28) 



The parameters resulting from the fit are listed in Table [2] 
Figure [8] shows the resulting M s - o~n/(n) relation. The Mach 
numbers were calculated as described earlier in Section 15.2. 11 
We show the relation only for the case when the cloud is defined 
with the Ay = 7 mag box, because the PDFs of the ellipsoidal 
definition can suffer from incompleteness (see Section 3.2 1. 

The advantage of using a log-normal model to estimate the 
shape parameters of the PDF is that allows estimation of the to- 
tal standard deviation of column densities, and from therein, di- 
rectly connect the shape parameters to the theory of turbulent 
density fluctuations (see Section 5.2 1. Under this scheme, we 
can also combine the yV( s - o~ni(n) measurements of the IRDCs 



with additional data points from K09 (listed in Table [3] calcu- 
lated from the data presented in Table l[^] of K09), who mea- 
sured the standard deviation of column densities in a sample 
of nearby (D < 200 pc) clouds. K09 data allow measurements 
of the standard deviation of the column densities in the range 
1 mag < Ay < 25 mag. Importantly, the standard deviation 
measured at this range measures reasonably accurately the stan- 
dard deviation measured from any wider range, assuming that 
the shape of the PDFs is even approximately correct. This is be- 
cause the probability of having column densities above Ay > 25 
mag is, in fact, so small that the contribution of those column 
densities to the standard deviation becomes negligible. For ex- 
ample, assuming that the true underlying PDF is similar to that 
of Ophiuchus (see Fig. ffl, but extends onwards from A v > 25 
mag in a power-law (or log-normal) fashion, the standard de- 
viation calculated from the range 1 mag < Ay < 25 mag is 
only 5% smaller than the standard deviation calculated from the 
all data above Ay > 1 mag. In summary, we can consider the 
K09 data compatible with the measurements of o"n/<n>> because 
the total standard deviation seems to be reasonably well recov- 
ered by them. However, we acknowledge that the measurement 
is based on different approach than the one we use in the case of 
IRDCs. It remains to be explored with a larger sample of IRDCs 
which exactly are the conditions under which the assumption of 
the compatibility is valid. 

We show in Fig. M the M s - o~n/(n) values for 14 nearby 
clouds with red triangles, calculated using the K09 PDF data and 
|Dame et al.| ( |200l) CO data. We list in Table[3]the Mach numbers 
derived for nearby clouds that span the range /W s ~ 7 - 22. The 
temperature of T = 10 K was adopted in the derivation. We also 
show in the plot the values derived by Brunt (2010]! and Padoan 
|et al.| ( fl997b| l that fall well into the range of our measurements. 
We note that the velocity dispersions for nearby clouds are de- 
rived from I2 CO data that generally probes slightly larger scales 
than 13 CO data. However the effect of this is relatively small; we 



3 We note that we have subtracted the contribution from noise to the 
standard deviations reported in K09, which typically is about 0.3 mag. 
We also note that K09 reports an erroneous for the value of Ophiuchus. 
Re-analyzing the same data, we get the standard deviation cr « 1.9 mag. 
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Table 3. Sonic Mach numbers derived for nearby clouds 



Cloud 


M s 


data 

N/(N) 


Ophiuchus 


10 


1.1 


Taurus 


10 


0.9 


Chal 


8.7 


1.0 


Chall 


12 


1.2 


CrA Cloud 


6.8 


1.1 


LDN1228 


11 


0.8 


LDN204 


9.6 


0.7 


LDN1333 


14.4 


0.6 


LDN1719 


8.4 


1.1 


Cha III 


11.5 


0.8 


Orion A 


21 


1.5 


Per Cloud 


11 


1.0 


Orion B 


17 


1.2 


Cepheus A 


22 


1.1 


California 


18 


0.7 



calculated that in Perseus, Ophiuchus and Taurus the mean dif- 
ference between the 12 CO and 13 CO linewidths is a factor of 1.2. 
We made an experiment in which we multiplied the linewidths 
of IRDCs by this factor and it does not significantly affect the 
relations derived in this section. 

Again, judging whether the M s - ct^kn) data of our IRDC 
sample are correlated or not suffers from the low number of 
clouds. If we treat the measurements for IRDCs and nearby 
clouds as equals (i.e., assume that both approaches trace reason- 
ably well the total column density fluctuations), the Pearson's 
correlation coefficients for the data are 0.61 and 0.66 for the 
Gaussian and fixed range velocity dispersions. This coefficient 
indicates correlation of variables with high significance (degrees 
of freedom = 20, p < 0.001). We also fitted a linear model to the 
data again using the linfitex procedure. In this calculation, we 
used the relative errors of 30% for the M s values, and 30% for 
the <tn/(n) values (the error of ctnkn) originates almost entirely 
from the uncertainty in the mean jj. We use the typical value 
given by the linfitex procedure). The slopes resulting from this 
were a i = 0.051 ±0.016 and a i = 0.042 + 0.013 for the Gaussian 
and fixed range velocity dispersions, respectively. Thus, the de- 
rived slopes are approximately 3-cr offset from zero. The error- 
weighted mean of the slopes is a\ = 0.047. We note that choos- 
ing the temperature of T — 15 K for nearby clouds (uniformly 
with the value chosen for IRDCs) instead of T - 10 K does 
not affect the significance of the correlation and results to simi- 
lar slopes (a; = 0.050 ± 0.015 and a x = 0.039 ± 0.012 for the 
Gaussian and fixed cases). 

We note that the uncertainty of our adopted Kg/i& ratio is ex- 
pected to play only relatively minor role in the uncertainty of the 
M s — (TM/(pf) relation. This is because we correlate specifically 
the mean-normalized column density variances with M s . Even 
though possible use of an erroneous opacity-law between NIR 
and MIR reflects directly to the column densities, their effect to 
mean-normalized column densities should be a second-order ef- 
fect. This is because, to first order, the column densities used in 
the PDFs (i.e., Ay > 7 mag) depend linearly on the opacity-law: 
for example, increase of k% by a factor of 1.5 would lower the 
maximum column densities approximately by that factor. And 
the PDF of a mean-normalized variable is invariant to multipli- 
cation of the variable by a constant factor. However, the uncer- 
tainty of the opacity-law does introduce a second-order effect to 
the shapes of the PDFs, because the response of column densities 
to the opacity-law is not exactly linear. Below Ay < 10-20 mag 
the large-scale component of the column density field becomes 



significant compared to the smaller-scale component and the sig- 
nal of the combined column density mapping technique (by con- 
struction, see Section 3.1 1 switches from being MIR-dominated 
to NIR-dominated. We illustrate the response of the combined 
column densities to the adopted k^/k^ ratio in Fig. C.l (online 
only). In this test, we kept the NIR-opacity constant and varied 
the MIR opacity by 15%, 30%, and 60%. The experiment shows 
that in these cases the PDF shape of even a mean-normalized 
column densities below Ay < 10-20 mag could be altered by 
the combination procedure. However, we include in the PDFs 
only column densities above Ay > 7 mag, and therefore it seems 
justified to assume that this kind of bias is clearly a second-order 
effect. For example, in the test above (which is for Cloud A), the 
<x N <n> changes about 5% and 20 % for the cases where MIR 
opacity is underestimated by 30% and 60%, respectively. 

Based on the experiments described above, we conclude that 
the observational data tentatively suggest correlation between 
the sonic Mach number and column density dispersion in molec- 
ular clouds. These measurements are, to our knowledge, the first 
direct measurement of the relationship of these variables and 
also the first direct detection of the correlation. However as also 
described above, we recognize the uncertainty of our measure- 
ments, rising mainly from the very small number of clouds in 
our sample and the uncertainties in defining "the cloud area" in 
our sample that is not tailored for the measurements of this par- 
ticular quantity. It would require a larger sample of IRDCs to 
strengthen this conclusion. In practice, such a study would need 
to target clouds in locations that are somewhat offset from the 
Galactic plane and as isolated as possible. 

Finally, we briefly speculate about the interpretation of the 
slope value we derive from the Ai s - ct^kn) data. As noted 
earlier, [Brunt et al.| ( |2010[ ) examined the 3D-2D variance ra- 
tio in a sample of isothermal turbulence simulations and mea- 
sured the ratios to be between R = [0.03,0.15]. If we simply 
assume that the ratio is in this range, the mean slope gives the 
values b = 0.047/ V{0.03,0.15) = {0.12,0.27} with the mean at 
b = 0.20. We report the 3-cr confidence range for this mean as 
the 3-cr deviation from the highest and lowest possible values: 
upper limit is given by b = (0.051 + 3 x 0.016)/ VO03 = 0.57 
and the lower limit by b = (0.042 - 3 x 0.016)/ V015 = -0.02. 
Thus, we report as our measurement b 



0.20 + ° g. This result 



points toward a somew hat lower value of b th an derived by Brunt 
( 120 10| (b = 0.493 and |Padoanetai] ( |l997b| l (b ~ 0.5). Low val- 
ues of b are generally predicted by (non-magnetized) isothermal 
simulations in which the turbulence driving is purely solenoidal 
(b = 1/3, e.g., |Federrath et al.|2008bl[20T0) l, while increasing the 
fraction of compressive modes increases the value of b toward 
unity (e.g., Federrath et al. 20 1 0| > . The presence of magnetic 
fields also modifies the M s - o~p/(p) correlation, and in general, 
marginally flattens the slope (lowers the b) compared to non- 
magnetized cases (e.g., Price et aT]|201 1) . Also, strong density 
dependence of the magnetic field makes the density fluctuations 
more weak l y dependent on M s (e.g., Ostriker et al.|2001 



et al.|20lT] Molina et al.|2012} . Even further, when the density 



Price 



dependence of the magnetic field strength approaches Sap, 



4 Note that [BTunt|j2010f reports the value b = 0.48, but they use an 
erroneous Mach number of 17 for Taurus. The c orrect Mach number for 
Taurus based on the data they use is A4 S = 10. [B7unt|p(jr0l l also uses a 
correction factor to account for column density variance not included in 
their measurement. If we do not use the correction factor, in which case 
their measurement of ctnkn) is very similar to ours, the value <xn/<n> = 
0.84 follows. Together with their measurement of R = 0.029 this leads 
to the value b = 0.49. 
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the M s - cthkn) 1S predicted to show no correla tion at Mach 
numbers greater than M s > 7 ( Molina et al.|2012 l. In this con- 
text, the fact that we detect a M s — <7"n/<n> correlation up to high 
Mach numbers supports the picture in which the magnetic field 
depends relatively weakly on density. This has also been sug- 
gested throu^hZe^rrmneffect observations (B oc p - 65 at n > 300 



cm 



-3 



Crutcheretal.|2010 i 



5.3. Comparison with earlier work 

How do our results regarding the M s - cr p /( p ) relation compare 
with the previous works on the topic? The Ai s - o>/(p> relation 
(or A1 s - o~nkn) relation) has not been widely examined by ob- 
servations. The most comprehensive work to our knowledge is 
that of Brunt (2010 ) who exploit CO line emission data and near- 
infrared dust extinction data to estimate M s and o-^kn), respec- 
tively, and use the method of Brunt et al. (2010 ) to infer cr p ^) 
and the proportionality constant b. Their observations indicate 
the value of b — 0.49 (see Footnote A similar technique was 
applied earlier to the dust extinction data of IC5146 by Padoan 
|et al.| ( [T997b| ) who estimated the value b * 0.5. Both of these 
works, however, only measure one M s — <Tnkn) ratio pair, and 
not explicitly the correlation between the two variables. 

|Goodman et al.| ( |2009| ) examined the PDFs and variances of 
the column density field in six sub-regions of the Perseus molec- 
ular cloud. The data (and the method) used by their work were 
very similar to what we employ in this paper. They tested if M s 
increases with increasing 2D variance (and mean) and they con- 
clude that it does not. Thus, the result of Goodm an et al.| f2009) 
suggests no correlation between o- p /(p) and M s . In the scope of 
the data we show in Fig. |7) it seem s clear that the contradiction 
between our and Goodman et al. (2009 ) results originate from 
the fact that the velocity dispersions in the sub-regions within 
Perseus span only a narrow range (Goodman et al. do not quote 
the linewidths of each sub-region, but the integrated linewidth 
over the entire Perseus is » 1 km s _1 , inevitably restricting the 
values to a rather narrow range). This leaves little leverage for 
the correlation to be detected, given the large scatter in the data 
shown in Fig. [8] 



6. Conclusions 

We present in this paper a new technique for probing the col- 
umn density structure of IRDCs. The technique is based on com- 
bining column density maps derived from NIR observations of 
background stars shining through the IRDCs (UKIDSS/Galactic 
Plane Survey data) and from MIR (8 pm) shadowing features 
against the Galactic background (Spitzer/GLIMPSE survey). 
We demonstrate the use of the method by deriving column den- 
sity data for ten IRDC complexes. We then use the resulting 
maps to examine the relation between the sonic Mach number 
and column density variance in the clouds. The main conclu- 
sions of our work are as follows. 

1 . The technique we present provides temperature-independent 
column density data over a unique combination of spatial 
resolution (~2") and dynamic range (Ay * 1 - 100 mag). 
Unlike the mid-infrared technique alone, the new method is 
well-calibrated at column densities (1 mag < Ay < 10 mag). 
The new information on the low-column density surround- 
ings of the IRDCs can significantly improve our understand- 
ing of the impact of the envelope material to the physics of 
the IRDCs (e.g. |Hernandez et al.|2012[ ). The technique is ap- 



20 +0 - 37 



plicable to clouds up to about D < 8 kpc, thereby potentially 
reaching a great number of IRDCs. 

The column density PDFs of the IRDCs in the range Ay = 
7-100 mag are, on average, consistent with a wide log- 
normal shape that has a width of cri nN ss 0.9. This shape 
is quite similar, although slightly shallower, to the PDFs of 
nearby active star-forming clouds reported by K09 in the col- 
umn density range in which the data overlap, i.e., Ay = 7-25 
mag. 

The IRDCs contain relatively high amount of high-column 
density material, comparable to the fraction observed in 
nearby active star-forming clouds such as Orion A, and pos- 
sibly even exceeding it. The average cumulative distribution 
function of IRDCs between 7 mag < Ay < 100 mag be- 
haves exponentially, CDF oc e -o.oixN j n me f ramewor k of 
the recent work on the star-forming rates of molecular clouds 
( |Lada et aL]|2012[ ), this suggests that the star-forming rates 
(per unit cloud mass) of IRDCs are at least comparable to 
active nearby clouds. 

We present the first direct observational measurement of the 
relation between the turbulent energy (sonic Mach number, 
Ai s ), and density fluctuations (column density dispersion, 
o~n/(n)), in molecular clouds. We report a tentative detec- 
tion of correlation with about 3-cr confidence. Linear model 
fits to the (M s ,(Tff/(N)) data yields the slopes of ctnkn) - 
({0.051,0.042} ± {0.016,0.012}) x M s for the two methods 
we use for measuring A1 s . 
Our results suggest the correlation coefficient b 
between the sonic Mach number and volume density vari 
ance, cr p i(p), in molecular clouds. In this expression, the 
quoted uncertainties indicate the 3-cr confidence interval as- 
suming that the 2D-3D variance ratio R is between R - 
[0.03,0.15] (as measured by |Brunt et aT1[20T0| . The value 
of b we derive is in agreement with earlier derivations by 
|Padoan et ah] ( |1997b| ) and |Brunt| ( |2010| ), that are, however, 
based on measurements of only single (A4 s ,o~n/(n)) ratio 
pairs. 

When combined with recent numerical and analytical predic- 
tions (Price et al.|201 1[ [Molina et al.|2012| ), the detection of 
correlation between Ai s and <tn/(n) ls suggestive of relatively 
weak density dependence of the magnetic field strength in 
the clouds. This is in agreement with the field strength mea- 
surements of Crutch er et al.| ( |20T0) ). 



While the conclusions given above remain tentative, the 
technique we present provides a good basis for measuring the 
relationship between M s - ctn/(n) in a statistically significant 
sample of molecular clouds in the future. 
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Fig. 8. Standard deviation of the mean-normalized column densities as a function of the sonic Mach number in eight IRDCs (black 
diamonds) and 14 nearby molecular clouds (red triangles). To compute the <tn/(n) values of the IRDCs, we fitted log-normal 
functions to the PDFs of the clouds and used the mean column density resulting from the fit to normalize the column density values. 
The PDF was constructed from the box-like area around the cloud as explained in Section 4.2 The 0-^/^ values of the nearby 
clouds are calculated over a dynamic range that should represent well the total column density dispersion (see text). The value 
derived by Brunt| (2010i and Padoan et al. ( 1997b| l is shown with a red asterisk. The dotted lines show linear fits to all data points. 
Left: yV( s is calculated from Gaussian fits to the mean spectra. Right: /V( s is calculated from a fixed velocity range (see Section 2.4 1. 
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Fig.A.l. The mean '3CO spectra of the IRDCs of our sample. The spectra have been averaged over the Ay > 7 mag box (see 
Section 3.2 for the detailed definition). The dotted vertical lines indicate the velocity interval chosen to represent the cloud. The 
red line shows a fit of a Gaussian to this velocity interval. The blue line shows another Gaussian fit, performed over the intervpj 
Vpeak - 1-5<T, Vpeak = 1.5ct where cr is the dispersion from the first Gaussian fit. The dispersions are shown in the panels. The third 
dispersion value, cr d gives the standard deviation of the data within the chosen velocity interval. 
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Fig. B.l. High-dynamic-range column density map of the cloud A, derived using a combination of MIR and NIR data. The green 
circle shows the largest ellipse from the Simon et al. (2006 ) catalog in the region. The white box outlines the region which was used 
alongside with the A v = 7 mag contour to define the region that is included in the analyses presented in Section[4] The contours are 
drawn at Ay = [7, 40] mag. The scale bar shows the physical scale assuming the distance as given in Table[TJ 
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but for the cloud D. The rectangular empty area results from missing NIR data. 
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Fig. B.6. Same as Fig. B.l but for the cloud F. The areas of missing data (marked with zeros) result from strong MIR nebulosity 



that hinders the MIR mapping technique. 
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Fig. B.8. Same as Fig. B. 1 but for the cloud H. 
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Fig. B.9. Same as Fig. B.l but for the cloud I. The areas of missing data (marked with zeros) result from strong MIR nebulosity 
that hinders the MIR mapping technique. 
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Fig. B. 10. Same as Fig. B.l but for the cloud J. The areas of missing data (marked with zeros) result from strong MIR nebulosity 
that hinders the MIR mapping technique. 
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Fig. C.l. Response of the column densities derived using the combined NIR+MIR technique to the adopted NIR-to-MIR opacity- 
law. The curves show the ratio of column densities derived with alternative opacity-laws (k^/k^ — nx 0.29) to the one adopted in 
this paper (k&/k^ = 0.29). The curves were calculated by keeping the NIR opacity fixed and by changing the relative MIR opacity. 
The solid blue, green, and red curves show the case in which the true opacity-law is 15%, 30%, and 60% higher than the adopted 
value (ks/Kk = 0.29). The dashed curves show the corresponding curves for the cases in which the ratio is lower. At Ay < 10 mag 
the technique relies dominantly on the NIR data, and therefore the changes in MIR opacity do not greatly affect the column density. 
The large-scale background component that is filtered out by the MIR data (and assumed to be recovered by the NIR data) has its 
maximum at Ay = 10-20 mag. Thus, at that range there is a transition from NIR-dominated to MIR-dominated regime. Finally at 
Ay > 20 mag the large-scale column density component starts to be small compared to the total column density, and thus the ratio 
approaches the value by which the original opacity law was modified (i.e., factors 1.15, 1.3, and 1.6). At extinctions lower than the 
10 mag threshold value, there is a small bias in the extinctions. It is caused by the fact that the background correction value (see 
Section 3.1 1 for some pixels can result from interpolation from neighboring pixels instead of from the difference of the NIR and 
MIR extinctions. The pixels for which this interpolation is performed are not necessarily the same in the cases where n — 1 and 
n + 1. 
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